Retardation of Particle Evaporation from Excited Nuclear Systems 

Due to Thermal Expansion 
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ABSTRACT 

Particle evaporation rates from excited nuclear systems at equilibrium 
matter density are studied within the Harmonic-Interaction Fermi Gas Model 
(HIFGM) combined with Weisskopf 's detailed balance approach. It is found 
that thermal expansion of a hot nucleus, as described quantitatively by 
HIFGM, leads to a significant retardation of particle emission, greatly ex- 
tending the validity of Weisskopf's approach. The decay of such highly 
excited nuclei is strongly influenced by surface instabilities. 
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I. INTRODUCTION 

Studies of nuclear heavy-ion reactions at Fermi bombarding energies (a few tens of 
MeV per nucleon) have faced conceptual and experimental challenges. This is the domain 
in which the low-energy dissipative reaction mechanism |l| is expected to morph into 
one that is characterized by a merging of time scales for collective motion and relaxation 
processes and the population of a new phenomenological realm |2J. For example, for 
nuclear temperatures of T ~ 5MeV, nucleon evaporation times of the order of t evap » 
20 fm/c have been estimated . Such short evaporation times describe a "prompt" 
nuclear decay occurring during the nuclear interaction between projectile and target, 
thus preventing a complete equilibration of the system. Furthermore, mechanical and 
chemical instabilities 4| expected for highly excited nuclear systems should exceed the 
range of conventional statistical models and cause novel nuclear decay modes ^,0]. 

On the other hand, experimental observations of the particle emission patterns 

in reactions at Fermi energies, induced by heavy projectiles, demonstrate the persistence 
of sequential decay of hot projectile-like fragments and their target-like reaction part- 
ners. Apparently, the primary hot fragments emit nucleons, light particles, and complex 
nuclear clusters with patterns characteristic of thermalized sources. While the presence 
of an additional, "intermediate-velocity source" component in the particle emission pat- 
terns (most likely, of dynamical origin) ^,11, 3| tends to complicate the analysis, the 



existence of intense sequential evaporation components demonstrates unexpectedly high 
(meta-) stability of the primary hot projectile and target-like nuclei. 

The present work investigates conditions for a consistency of the apparently long 
particle emission times with high nuclear excitations. It is shown that, consistent with 
the above observations, highly excited nuclear systems are indeed similar to a classical 
compound nucleus, once their expansion to the equilibrium matter density is taken 
into account. The particle emission timescales are then long enough to allow for a 
meaningful equilibration of the system in the sense of Weisskopf 's detailed balance jl^ . 
Relative decay rates for various statistical decay channels are governed by the entropy at 
a respective "transition" state, in which the state of the surface plays an important role. 
The transition state may be identified with an emission barrier for particle emission and 
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a separation saddle for statistical fission-like decay. 

In the following section, the Weisskopf evaporation model is briefly reviewed. The 
effects of nuclear expansion are described in Section III, followed by a presentation of 
calculations in Sect. IV and brief conclusions (Sect. V). 



II. WEISSKOPF'S PARTICLE EVAPORATION TIME-SCALES 

The evaporation time scales for (metastable) excited nuclear systems can be con- 
veniently evaluated using Weisskopf's detailed balance [yj approach. This approach 
relates the probability W(E A , e) for the emission of a particle with energy e by a com- 
pound nucleus with mass number A and excitation energy E* A to the cross section cr(e) 
for the reverse process of absorption of such a particle by a nucleus of mass number B 
and corresponding excitation energy Eg, 

TO«)*=»(«e^ (i) 

7T ft U A {Ea) 

where uja and ujb are the level densities of nuclei A and B at respective excitation 
energies, m is the mass of the emitted particle, ft is Planck's constant divided by 2tt, 
and g denotes the spin degeneracy. 

In Eq. [TJ excitation energies E* A and Eg of the compound nucleus A and the residue 
B are related via energy conservation 

E *B = E A ~ e - Qgg, ( 2 ) 

where Q gg is the Q — value for the particle separation from the ground state of nucleus 
A, leaving the residue B also in its ground state. 

By expressing level densities uo in terms of entropy S 



u A {E A ) = e SA{E ^ and u B (E* B ) = e Ss ^) (3) 
one rewrites Eq. Q in a form 



ir 2 h 



(4) 
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which is well suited for use in conjunction with models such as HIFGM (k^ 15]. Ulti- 
mately, the decay width T for a particular particle channel is obtained by integrating 
Eq. H]over particle energy e and then by multiplying the result by H, i.e., 

T(E* A ) = h r W{E% e)de = OIL [°° a (e)ee s ^- Q ^- SA ^de. (5) 
Jo n 2 n Jo 

To be able to make use of Eq. EJone needs additionally a prescription for calculating 
the cross section as a function of bombarding energy. In the context of Weisskopf's 
approach, one assumes customarily that the absorption cross section is equal to the 
geometrical cross section of the nucleus. This corresponds to a constant value of a = 
Txr 2 for neutrons and an energy- dependent, Coulomb corrected expression for charged 
particles, 

a = vrr 2 (l-^), (6) 

where r Q is the nuclear radius parameter and Vc is the magnitude of the Coulomb barrier. 
Obviously, Eq. El is suited for both, charged particles and neutrons, with Vc — for the 
latter. 

A useful approximation for the right-hand side of Eq. 03 can be obtained ^| for the 
neutron decay width. It assumes an identical functional dependence of the entropy on 
thermal excitation energy for both, mother and daughter nucleus, i.e., S a (x)=Sb(x). 
Then, expanding the entropy Sb in a Taylor series about E* A to first order results in 

Sb(E* a - Q gg - e) - S A {E* A ) = * (Q 99 + e) = -^(Q gg + e), (7) 

dE A 1 A 

where use was made of the fact that the derivative of the entropy with respect to ex- 
citation energy, dS A /dE A , is equal to the inverse of the temperature, T A , of nucleus A. 
Thus, for the neutron decay width one obtains an approximate expression 

UE* A ) = ^r D 2 Tie-%\ (8) 
Tin 

revealing the strong temperature dependence of the decay width. It is this exponential 
temperature dependence that has raised concerns regarding the applicability of Weis- 
skopf's equilibrium approach at higher excitation energies. Rationale given in Sect. Ill 
will demonstrate a much weaker effective temperature dependence of the decay width. 
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In figures shown further below, the decay rate is expressed in terms of the decay time 
scales parameter r, which is related to the decay width T via Heisenberg's uncertainty 
principle, r = fi/Y. 



ILL EFFECTS OF THERMAL EXPANSION OF NUCLEAR MATTER 

Like most liquids, nuclear matter is expected to respond to thermal excitation by 
expanding, so as to arrive at a state of maximum entropy for any given total excitation 
energy. While thermal expansion should be predicted by any realistic model of bound 
nuclear matter, a particularly simple expression for the dependence of the asymptotic, 



jv the Harmonic In- 
13 . Il^ l are described 



equilibrium matter density on total excitation energy is provided 
teraction Fermi Gas Model (HIFGM). The essentials of this model 
below. 

The central notion of the HIFGM is entropy. The relation between this entropy and 
the thermal excitation energy El herm is assumed to be described by the regular Fermi-gas 
expression, such that 



S(E; ot ) = 2^a(E? ot -E* compr ), (9) 

where E^ ot is the total excitation energy and E* ompr is its collective compressional (poten- 
tial energy) part. The difference between these latter two energies defines the random, 
thermal excitation E^ herm and, hence, the entropy. Furthermore, the HIFGM adopts the 
matter density dependence of the level density parameter a ("little-a") germane to the 
Fermi-gas model, 

a = a (— )-f, (10) 

Po 

where a Q is the ground-state value of the level-density parameter and p Q is the ground- 
state matter density. For simplicity, the HIFGM assumes a quadratic dependence of the 
compressional energy, E* ompr , on the relative matter density [lfil ] 

Ecompr - E bind (l - "J") 2 \ (11) 
Po 

where Eb in( [ is the ground-state binding energy of the system. 

Given Eqs. I9T ITUI the equation for maximum entropy as a function of relative density 
can be resolved analytically, yielding jl^ 
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— = t( 1 + \/ 9 -# L )- ( 12 ) 

Assessing qualitatively the effects of thermal expansion on the decay width (see Eq. |HJ), 
one notes that, according to Eq. El the matter density decreases with increasing total 
excitation energy. Consequently, the compressional part of the excitation energy and the 
level density parameter, both increase. The latter two trends reduce the rate of increase 
of the system temperature with increasing excitation energy, as compared to that without 
expansion. This follows from the Fermi-gas model relationship between temperature 
(T), level-density parameter (a), and thermal part (E* herm ) of the excitation energy, 

E therm = ^tot ~ ^ compr = a ^ ■ (13) 

This effect can be viewed as a manifestation of the LeC hotelier Principle. The reduc- 
tion in effective temperature from the value that it would have at normal matter density 
leads to a reduction of the decay width, as approximated by Eq. |HJ This reduction is 
partially offset by the effects of a) an increased value of the inverse capture cross section 
(cr) and b) a neutron binding energy (Q gg ) that is reduced by an amount equal to the 
compressional energy per nucleon. Results of calculations performed using the "exact" 
Eq. are discussed in the next section. 



IV. RESULTS OF CALCULATIONS 



Results of model calculations for the Weisskopf particle evaporation time scales are 
shown in Figs. 1 and 2. In Fig. 1, the evolution of the average nucleon evaporation 
time r with total excitation energy per nucleon, predicted for expanded Fermi matter 
at equilibrium density (solid line), is compared to that resulting for Fermi matter forced 
to stay at ground-state density (dashed line). As seen in this figure, if allowed to 
expand to equilibrium density, a nucleus excited to E^ ot /A> 8 MeV regains a degree 
of (meta)stability against nucleon evaporation comparable to that of a ground-state 
density nucleus at the much lower excitation of E^ ot /A^ 4.5 MeV . The small increase 
in r predicted to occur near the high-excitation end of the curve is related to the negative 
heat capacity predicted by the HIFGM in that energy domain. Both these effects are 
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FIG. 1: Excitation-energy dependence of Weisskopf 's nucleon evaporation time scales for nuclei 
at equilibrium (solid line) and ground-state (dashed line) matter density. The dotted line 
illustrates results of calculations with T = absorption cross section. 

due to the fact that the considered maximum in entropy is not absolute, but conditional, 
obtained by imposing a certain functional form on the radial distribution of the nuclear 
matter density. 

In Fig. 2, partial time scales are illustrated separately for neutron emission and for 
proton emission. In both cases, thermal expansion is seen to lead to a retardation of 
emission, extending the range of applicability of Weisskopf 's approach. 



8 
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FIG. 2: Excitation-energy dependence of Weisskopf 's neutron (bottom curves) and proton (top 
curves) partial evaporation time scales for nuclei at equilibrium (solid line) and ground-state 
(dashed line) matter density. 

V. CONCLUSIONS 

The present calculations have shown that thermal expansion of a nucleus to equilib- 
rium density has a stabilizing effect on the excited system, such that the most likely 
statistical decay channels are suppressed. This may give the system time to reach a 
more complete statistical equilibrium, where various surface modes are excited as well. 
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Such surface modes may be viewed as doorway states associated with cluster emission. 
It has been demonstrated elsewhere that, when equilibrated at high excitation ener- 
gies, such fission-like decay modes may compete successfully with nucleon emission. It 
should also be kept in mind that, with every neutron that escapes before these cluster 
emission modes are equilibrated, the probability for the emission of subsequent neutrons 
decreases, due to isospin effects, until it equals that for proton emission. This isospin ef- 
fect adds additional credibility to the idea that Weisskopf 's approach may be applicable 
at much higher nuclear excitation energies than commonly thought. 

It is worth noting that, in its present schematic form, the HIFGM assumes that 
thermal expansion occurs in a self-similar fashion, i.e., that this expansion can be reduced 
to a rescaling of the radial coordinate. In addition, for the sake of simplicity, it is assumed 
that the effective nucleonic mass is equal to the free nucleon mass. A more complete 
modelling 17, Q] calls for more thorough and detailed treatment of the diffuse nuclear 



surface and its evolution with increasing excitation. In such a treatment, it would be 
highly desirable to account for the finite range of the effective nuclear interaction. This 
finite range may alter the thermostatic properties of surface matter with respect to those 
of bulk matter at equal density and temperature. Furthermore, it is necessary to include 
in a comprehensive modelling also thermal surface oscillations of the hot nuclear system, 
as the role of such oscillations is expected to increases [13] with increasing excitation 
energy (through decreasing surface tension or increasing surface entropy). It is to be 
expected that, in line with LeC hotelier' s Principle, both, the bulk nuclear matter and 
the surface region will evolve in a concerted manner that increases entropy and stabilizes 
the overall system with respect to single nucleon emission even further. 
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